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The velocity-dependent cusp anomalous dimension is by definition the ultraviolet (UV) anomalous dimension 
of a Wilson loop with a cusp. It appears in many physically interesting processes. In this talk we present recent 
progress in the analytic calculation of that quantity. 



1. INTRODUCTION 



Wilson loops are very fundamental and impor- 
tant quantities in gauge theories [1 . They are es- 
sential for defining nonlocal gauge invariant quan- 
tities. They also contain information about lo- 
cal operators via the operator product expansion. 
Moreover, it turns out that certain Wilson loops 
defined for specific contours appear as an effective 
description of certain physical processes. 

The velocity-dependent cusp anomalous di- 
mension is a case in point. It is defined as the ul- 
traviolet (UV) anomalous dimension of a Wilson 
loop with a cusp. It appears in many physically 
interesting situations. For example, in heavy 
quark effective theory (HQET), it describes the 
infrared (IR) divergences of massive form factors 
and scattering amplitudes, see e.g. |2)3|4|5j . 

In a M = 4 super Yang-Mills (SYM), it is nat- 
ural to define the locally supersymmetric Wilson 
loop operator |6|7j 



Pexp 



\dx\n • $ 



(i) 



where n is a vector on S 5 . It parametrizes the 
coupling of the Wilson loop to the six scalars <S 
of the theory. We consider as the integration con- 
tour a cusp formed by two segments along direc- 
tions (momenta) and and allow the two 
segments to couple to the scalars through vectors 
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Figure 1. A Wilson line that makes a turn by an 
angle <p in Euclidean space. The two segments are 
directed along p M and respectively. The vec- 
tors ft and n' are internal vectors that determine 
the coupling to the scalars <l>, see eq. |l]). 



n and n', see Fig. [Tj Then, the vacuum expecta- 
tion value (W) of the Wilson loop will depend on 
the angles 



P'Q 



(2) 



as well as on the 't Hooft coupling A = g 2 N, and 
the rank of the gauge group N. 

If Auv an d Air are short and large distance 
cutoffs, respectively, then the divergent part of 
the vacuum expectation value of the Wilson loop 
takes the form |l|8j 



(W) ~ exp 



A uv r 

log - — r, 



A 



IR 



cusp 



(3) 



cusp anomalous dimension 



This defines the 

r cusp (0,M,iv). 

We will begin by describing the dependence of 
r cU sp on the angles (j) and 6 in section [2j and then 
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introduce, in section [3J the set of integral func- 
tions needed to express loop-level results. We will 
then review perturbative results to two loops, and 
present the new three- loop result. We present the 
light-like limit as a check in sectionJEJ and de- 
scribe our method of calculation in |6l We then 
present an exact result for r cusp at small angles 
in section [7J In section |8j we discuss a new scal- 
ing limit. Finally, we conclude and present future 
directions in section [9l 

2. KINEMATICS AND LIMITS 

For reasons that will become apparent 
presently, we will mostly be interested in the <ft 
dependence of the Wilson loop. It is convenient 
to introduce a new variable x = e 1 ^. The com- 
putation we are considering is invariant under 
(j) — )> —cj). This corresponds to an inversion sym- 
metry in x. 

Note that the dependence of r cusp on is sim- 
ple. It can only occur through Wick contractions 
of scalars, and because of SO (6) invariance it ap- 
pears only through ft • n' = cos 6. Therefore, at 
L loops, r cusp is a polynomial in cos#, of max- 
imal degree L. Having made this observation, 
we find that it convenient to introduce the vari- 
able £ = (cos 6 — cos <p) j (i sin <p) . Two important 
cases are 6 = (constant coupling to scalars) and 
= 7r/2 (scalars on opposite edges are orthogonal 
to each other), which lead to £ = (1 — x)/(l + x) 
and £ = (1 + x 2 )/(l — x 2 ), respectively. 

There are several special cases of the angles 
that are of particular interest. When the geomet- 
ric angle <p and internal angle satisfy (j) = ±0, 
the anomalous dimension vanishes. Note that £ 
vanishes in this case. For = 0, this corresponds 
to (j) — > 0, i.e. x —> 1, the case of a straight 
line. The small angle limit is related to energy 
loss of an accelerated quark, and is known ex- 
actly |10lllj . The limit (j) -> 7r, i.e. x -> — 1, 
is related to the quark- ant iquark potential. This 
limit is subtle and requires a resummation of cer- 
tain diagrams, see |12|13|14fT5] . 

The above limits can be defined in Euclidean 
space. There is an intrinsically Minkowskian limit 
that is also of interest. When <p — iip, <p ^> 1, 
i.e. x — >• 0, the cusp anomalous dimension di- 







x = e i( t > 








-1 
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Figure 2. r cusp is real for < x < 1, and be- 
low threshold, where x = e l€ ^ is a phase. Above 
threshold, we have — 1 < x < 0, with x having 
an infinitesimal imaginary part. The zigzag line 
denotes a branch cut along the negative real axis. 
Recall that r cusp has is an inversion symmetry 
x — ^ 1/x. 



verges linearly in the (/?, to all orders in the cou- 
pling constant [9]. The coefficient of the lin- 
ear divergence is the well-studied light-light cusp 
anomalous dimension; the latter can also be ob- 
tained from the anomalous dimension of high spin 
operators |16|17|18] . The Wilson loop approach 
considered here is a very efficient way of comput- 
ing this quantity. 

3. INTEGRAL FUNCTIONS 

Let us discuss the different kinematical regions 
for complex x. It is useful to recall the relation- 
ship of r cusp to IR divergences of massive form 
factors, such as 7* — > e + e~, which have the same 
analytical structure. 

r cusp is real in the Euclidean region < x < 1, 
and in the physical region below threshold (of cre- 
ating two massive particles), which is the unit cir- 
cle (i.e. for (j) real). The threshold is at x = — 1, 
and the physical region above threshold goes from 
— 1 to along the negative real axis. See Fig. [2] 
We find that all results can be expressed in terms 
of harmonic polylogarithms (HPL) [19]. These 
functions are generalizations of classical polylog- 
arithms, and appear naturally in this problem. 
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They are defined iteratively by 
H au a 2 ^ an {x)= [ f ai (t)H a2i ... tan (t)dt, (4) 

JO 

where the integration kernels are defined 
as fi(x) = (1 - x)" 1 , f (x) = x~ 1 , and 
f-i(x) = (1 -fx) -1 . The degree (or weight) one 
functions needed to start the recursion are defined 
as 

Hi(x) = -log(l-x), H (x) = log(x) , 

ff_i(aO = log(l + a?). (5) 

The subscript of H is called the weight vector. 
A common abbreviation is to replace occurrences 
of m zeros to the left of ±1 by ±(m + 1). For 
example, i^ ,o,i,o,-i(^) = #3,-2(2;). 

HPLs have simple properties under certain ar- 
gument transformations, and one can use their 
algebraic properties in order to make their asymp- 
totic behavior manifest. We refer the interested 
reader to ref. [19 j. A very useful computer alge- 
braic implementation has been given in ref. [20] . 

The perturbative results given in section [4] can 
be straight forwardly evaluated numerically in the 
region < x < 1. Other regions can be reached 
by analytical continuation, respecting the above 
branch cut properties. A related discussion is 
given in ref. [21 for the process gg — >• h via a 
massive quark loop. 

4. RESULTS UP TO THREE LOOPS 

The cusp anomalous dimension of the standard 
bosonic Wilson loop operator was computed in 
QCD at two loops in the pioneering paper [5j. 
This result was later simplified [22], and recom- 
puted in J\f = 4 SYM for the supersymmetric 
loop operator defined in eq. ([!]) in refs. |23)24j . 
The answer in Af = 4 SYM can be written as an 
expansion in the 't Hooft coupling A, 

Tcusp = ^2 (g^2 ) r cusp • ( 6 ) 

To two loops it is given by, 

r« p = -£iog*, (7) 

Ifulp = 1^ N X (log 2 X + 7T 2 )] 



+\e [#1,1,1 + 2jfi, 2 ] • (8) 

Here and in the following, the HPLs are under- 
stood to have argument 1 — x 2 . In ref. [14], we 

(3) 

computed the three-loop value Tcusp- Before out- 
lining the method of calculation in section [6] let 
us present the result. Here we write it in a very 
compact form, due to [25] , 

Ifulp = -^log^lGg^ + TT 2 ) 2 

-^ 2 [3C3#i,i + C 2 (2#i,2 + 4#2,i + 3ffi,i,i) 
3 11 

+ 2^1,1,1,2 + -#1,2,1,1 + 2#2, 1,1,1 + -^#l,l,l,l,l] 

+ 3 [4ff M>3 + 4# 1)2 , 2 + 4^1,1,1,2 

+2^1,1,2,1 + 2^,2,1,1 + 3#i,i,i,i,i] • (9) 

Let us now discuss these formulas. 

It is not hard to verify the expected inversion 
symmetry x —> 1/x, i.e. <\> — >> — 0, using standard 
relations between HPLs of related arguments [19] . 
Note that £ is antisymmetric under this transfor- 
mation. 

We see that at each loop order, we have func- 
tions of uniform degree of transcendentality (i.e. 
weight of harmonic poly logarithms) (2L — 1). 

It is remarkable that the result can be written 
using argument 1 — x 2 (or x 2 ) only. This is not the 
case in general for individual Feynman integrals. 
Moreover, we see that there are no relative signs 
within the coefficients of the different powers of £ 
in the expressions above. This property also holds 
at least up to L = 6 for the £ L terms at L loops, 
which have recently been computed analytically 
[25]. 

5. RELATION TO LIGHT-LIKE CUSP 
ANOMALOUS DIMENSION 

Let us use the above three-loop result to ana- 
lytically verify the relation 

lim r cusp = - l - log x r~ sp + O(x ) , (10) 

where T^ sp is the anomalous dimension of a light- 
like cusped Wilson loop [26 . The asymptotic 
limit x — >• is easy to perform on our results, 
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since the logarithmic behavior of HPLs can al- 
ways be made manifest [19]. From eqs. ([7]), ([8| 
and Q it is easy to verify its correct three-loop 
value [27] 



■ cusp 



87r 2 7 



0(A 4 ). 



(11) 



6. RELATION TO MASSIVE SCATTER- 
ING AMPLITUDES IN N=4 SYM 

The velocity-dependent cusp anomalous di- 
mension governs the infrared divergences of mas- 
sive scattering amplitudes. In M = 4 super Yang- 
Mills, one can introduce masses by giving a vac- 
uum expectation value to some of the scalar par- 
ticles. In this way, at the planar level, one can de- 
fine finite four-dimensional scattering amplitudes 
which have an exact dual conformal symmetry 
[28] . The simplest case is the four-particle scat- 
tering amplitude M±(u,v), which is a function 
of two conformally invariant variables u and v, 
which in turn are defined in terms of the momenta 
and masses. 

It was argued in |29I30I14| that the limit m«1 
of this function is determined by the cusp anoma- 
lous dimension in the following way, 



lim log M±(u, v) = log u T c 

u— )-0 



,(v) + O(u ). (12) 



The relation between v and the previously used 
variables is v = x/(l — x) 2 . 

Let us explain how we used this limit to com- 
pute r cusp at three loops [14 , starting from a 
known representation of the three-loop scattering 
amplitude |27|29] , We observed that all known 
form factor and Wilson line integrals of this type 
could be expressed in terms of harmonic polylog- 
arithms of argument x. Assuming that this holds 
for r cusp at three loops reduced the calculation 
to the problem of determining a number of coeffi- 
cients. We used Mellin- Barnes techniques [29] to 
compute the asymptotic limit x — ^ 0, keeping not 
only powers of logx, but also terms suppressed 
by powers of x. In this way, we were able to fix 
our ansatz completely. 
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Figure 3. The form factor (or Wilson loop) inte- 
grand is obtained by taking the limit where all 
dual integration points, x^x^x^ in (a), come 
close to the external point x<i- Here we show 
a three-loop integral for the four-particle ampli- 
tude (a) giving rise to contribution for the Wilson 
loop, shown in (b). Thick lines represent massive 
(or eikonal) propagators, and thin lines represent 
massless propagators. 



We note that recently, our assumption was 
proven for two infinite classes of integrals, and 
moreover that it is possible to directly evaluate 
the corresponding integrals [25] . 

We wish to emphasize that the relation to the 
four-particle amplitude discussed above can also 
be used at the level of the integrand, where the 
dominant region of integration as u —> corre- 
sponds to all dual loop moments^] approaching 
one of the external dual variables, say x%. In 
that region, one can simplify the integrand (to 
logarithmic accuracy) to obtain integrals of form 
factor type. In the same spirit, one can further 
approximate them in order to get a Wilson line 
representation. We illustrate this simplification 
in Fig. [3] 

In summary, this means that the integrand of 
the Wilson line calculation can straightforwardly 
be obtained at higher loop orders from the corre- 
sponding integrand of massless scattering ampli- 
tudes |31|32|33j , extended properly to the massive 
case [28 29 . This should represent a useful start- 



2 Here we are using the notion of dual (or region) variables 
for a planar graph. 
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Figure 4. Bremsstrahhmg function B(X) (solid 
black line) given in eq. (14). The dotted green 
lines on the left denote two-loop and three-loop 
approximations, respectively, and the dotted red 
line on the right corresponds to the first two terms 
of the strong coupling approximation. 



ing point for calculations at higher loop orders, as 
well as for finding and proving conjectures about 
the structure of r cusp , see e.g. section [8] 

7. EXACT RESULT FOR SMALL 
ANGLES 

It turns out that there is a special limit of r cuS p, 
namely the small angle limit \0 2 — (j) 2 \ <C 1, which 
can be computed exactly in (/>, A and N. Here we 
simply quote the result and refer the interested 
reader to [10 for more details. Specifying to the 
planar case N ^> 1, we have 



1 



1-0 2 /tt 



r£(A) + ..., (13) 



with A = A(l - <t> 2 /n 2 ), and 

i Vxh(Vx) , 



B(X) 



4^ 2 h(VX) 



G(l/N 2 



(14) 



where Ij is the modified Bessel function. 

Expanding this exact result to the third order 
in A, one reproduces the coefficient of £ in eqs. 
Q, (|8| and ([9|, respectively. 

One can also expand our exact result at strong 
coupling. This can be compared against taking 
the limit \(f) — 0\ <C 1 of the formulas in [24]. In 
both cases, we find 



(0-<p) 



27T ^ 



(15) 



For a constant coupling of the Wilson loop to 
scalars, i.e. = 0, eq. (TT3J) becomes 



■ cusp 



-0 2 5(A) + 0((/> 4 ), = 0. (16) 



This is related to the radiation emitted by a 
quark undergoing a small change in direction [10| . 
The "Bremsstrahlung" function B(X) is plotted in 
Fig.gJ 

8. SCALING LIMIT 

The existence of the parameter 6 allows us to 
define a new scaling limit. Recall that r cusp de- 
pends on cos in a polynomial way. It is easy to 
see that the terms with the highest power cos L 
at L loops come from a diagram with L scalar 
propagators ending on each Wilson line. We can 
isolate such terms by taking the limit [14] 



A H> , e t& oo , 



«= -r-zXe* fixed. (17) 

47T 2 V 7 

In this way, we get a non-trivial function of </>, 
k, and N. We see that at leading order (LO) in 
the limit, only ladder diagrams remain. We will 
denote the leading order approximation by r^ p . 

The sum of the ladder integrals satisfies a 
Bethe-Salpeter equation. Since we are only inter- 
ested in the leading UV divergence (of the sum of 
the ladder integrals), we can simplify the equation 
to a one-dimensional Schrodinger problem, 



1 



cosh y + cos ( 



*(y) = -^*(y).(18) 



Here y is a variable in an auxiliary space (related 
to relative positions on the two edges of the Wil- 



son lines), and T\ 



LO 

cusp 



-^o> where is com- 



puted for the ground state of the system. 
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In the case of zero angle, the potential in eq. 



([18]) becomes the integrable Poschl- Teller poten- 



tial cosh (y/2), and one can solve for 
closed form. One finds [14 



1- vTT^ 



16 



1 + vT" 



1 + ^ + 2^ 



(19) 



Moreover, for general angle, the problem can 
be solved perturbatively. One can show that to 
any order L, the LO answer for r cusp can be writ- 
ten in terms of harmonic polylogarithms of degree 
(2L — 1), and an algorithm to compute the so- 
lution was given [25]. Moreover, the result was 
given explicitly up to six loops. This calculation 
also confirmed analytically the £ 3 term at three 
loops that was given in equation ([9|. 

There are a number of interesting properties of 
these results [25]. A study of the explicit results 
to six loops shows that in fact only a subset of 
harmonic polylogarithms is needed. We already 
mentioned in section [4] the absence of relative 
signs. Finally, when studying the limit x — >• 
discussed in section [5j one finds, at least up to 
six loops, that the asymptotic expansion of T^® 
does not require multiple zeta values [34] of depth 
two or higher. 

One can also compare the solution of the 
Schrodinger equation to what one obtains from 
taking the scaling limit of the corresponding 
string theory result. A priori, the answer ob- 
tained does not have to be the same, since the 
order of limits could be important. However, the 
explicit calculation showed agreement [T4] . 

This analysis has recently been extended to the 
next-to-leading order (NLO) in the scaling limit 
[25] . (See also ref. [35 for an analysis of the 
special case <\> — >• tt.) One finds that at NLO there 
are two infinite classes of scalar integrals, shown 
in Fig. [5] They satisfy modified Bethe-Salpeter 
equations and can be computed algorithmically 




Figure 5. The two infinite classes of scalar inte- 
grals the contribute to r cusp at NLO in the scaling 



limit (17). The integral on the right contains a 
numerator factor (k± + /C2) 2 , where k\ and k^ are 
the momenta flowing along the propagators with 
the arrows, respectively. 



Mills. 

Recently, it was shown how to apply integra- 
bility techniques to this problem [36 37 38 . It 
would be interesting to reproduce the results ob- 
tained here using these methods. 

While the results presented here are valid for 
supersymmetric gauge theories, the structures 
found are also be relevant for QCD. Interesting 
extensions of our work would be the calculation 
of the three-loop value of the cusp anomalous di- 
mension in QCD, as well as the first non-planar 
correction, which appears for the first time at four 
loops. 
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9. DISCUSSION 



We have reviewed various new results for the 
cusp anomalous dimension in M = 4 super Yang- 
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